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Abstract: The decay constants of the pseudoscalar mesons B and are evaluated 
from QCD sum rules for the pseudoscalar two-point function. Recently calculated per- 
turbative three-loop QCD corrections are incorporated into the sum rule. An analysis in 
terms of the bottom quark pole mass turns out to be unreliable due to large higher order 
radiative corrections. On the contrary, in the MS scheme the higher order corrections are 
under good theoretical control and a reliable determination of Jb and Jb^ becomes feasible. 

Including variations of all input parameters within reasonable ranges, our final results 
for the pseudoscalar meson decay constants are = 210ibl9 MeV and fs^ = 244±21 MeV. 
Employing additional information on the product y/Bsfs from global fits to the unitarity 
triangle, we are in a position to also extract the S-meson 5-parameter Bb = 1-26 ± 0.45. 
Our results are well compatible with analogous determinations of the above quantities in 
lattice QCD. 
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1 Introduction 



Experimental efforts in recent years have provided us with a wealth of new information on 
decays of bottom hadrons. To achieve a good understanding of these data, also the impact 
of the strong interactions has to be controlled quantitatively. This requires the accurate 
calculation of hadronic matrix elements involving S-hadrons. Generally, hadronic ma- 
trix elements contain contributions from low energies and thus non-perturbative methods 
should be employed for their evaluation. Current approaches include lattice QCD, QCD 
sum rules and the effective theory of heavy quarks (HQET). In this work, we shall consider 
a calculation of the simplest type of hadronic matrix elements, namely the pseudoscalar 
B- and Bg-meson decay constants fs and fs^ in the framework of QCD sum rules 

The pseudoscalar decay constants parametrise 5-meson matrix elements of the axial- 
vector current with the corresponding quantum numbers and are defined by 

(0|(g7M75?>)(0)|S(p)) = tfsPf. and (0| (57^,756) (0)|5s(p)) = ^fB.P,. (1.1) 

Throughout this work we assume isospin symmetry and q can denote an up or down quark. 
Weak interactions induce the leptonic decay of the S-meson. For example Jb then appears 
in the decay width of the process bu —>■ Wi which takes the form 

T{B-^l-u,) = ^Ktl'fBmJmsfl-^], (1.2) 



completely analogous to the corresponding decay of the light pseudoscalar mesons. Despite 
the suppression by the small factors mf and iKitp, there is some hope that the leptonic 
decay B Wi can be measured at the B-factories within the next years. Once Jb is 
assumed to be known, this would provide a very clean determination of \Vub\- In any case, 
Jb is an important quantity for it also enters more complicated hadronic matrix elements 
of 5-mesons like form factors or matrix elements of four-quark operators. 

The calculation of heavy meson decay constants in QCD has a rather long history. 
For charmed mesons, they were first considered in 1^,^, whereas the extraction of J'b 
from QCD sum rules was investigated in |]7|-p3|. The first determination of fs [0] dates 



back already twenty years. Nevertheless, due to recent theoretical progress, we find it 
legitimate to reconsider this problem. Very recently, the perturbative three-loop order 
correction to the correlation function with one heavy and one massless quark has been 



calculated p6l,p!7| for the first time. It turns out that in the pole mass scheme, which 
was used for most previous analyses, due to renormalon problems |T^, the perturbative 
expansion is far from converging. However, taking the quark mass in the MS scheme [ff9| , 
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a very reasonable behaviour of the higher orders is obtained and a rehable determination 
of /b becomes feasible.^ 

The starting point for the sum rule analysis is the two-point function of two 

hadronic currents 

v^(p2) ^ ^ [dxe^P^{Q\T{Ux)j,{Oy}\Q), (1.3) 



where Q denotes the physical vacuum and j^ix) will be the divergence of the axialvector 
current, 

j5{x) = {M + m) ■.q{x)n5Q{x):, (1.4) 

with M and m being the masses of Q{x) and q{x). In the following, Q{x) denotes the 
heavy quark which later will be specified to be the bottom quark, whereas q{x) can be one 
of the light quarks up, down or strange. Note that the current ^5(0;) is a renormalisation 
invariant operator. In the case of {ub) the corresponding matrix element is given by 

(m6 + m„)(0|(m756)(0)|E) = /sml , (1.5) 

where ttlb is the S-meson mass. 

Up to a subtraction polynomial which depends on the large p"^ behaviour, \I^(p^) satisfies 
a dispersion relation (for the precise conditions see pi|): 

^(p2) = f rfs + subtractions , (1.6) 

J [s — p^ — ^0) 


where p(s) is defined to be the spectral function p(s) = Im\l/(s + i0)/7r. To suppress 
contributions in the dispersion integral coming from higher excited states, it is further 
convenient to apply a Borel (inverse Laplace) transformation to eq. ( |1.6| ) which leads to0 

00 

uBu^iji^) = u^{u) = J e-'/^pis)ds. (1.7) 



Bu is the Borel operator and the subtraction polynomial has been removed by the Borel 
transformation. As we shall discuss in detail below, the left-hand side of this equation is 
calculable in renormalisation group improved perturbation theory in the framework of the 
operator product expansion, if the Borel parameter u can be chosen sufficiently large. 



^After completion of this work, we became aware of an independent analysis on the same subject [ pO[ 
where also the new order corrections are included, however employing the framework of HQET. 
^AU relevant formulae for the Borel transformation are collected in Appendix A. 
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Under the crucial assumption of quark-hadron duality, the right-hand side of eq. ( |1.7] ) 
can be evaluated in a hadron-based picture, still maintaining the equality, and thereby 
relating hadronic quantities like masses and decay widths to the fundamental Standard 
Model parameters. Generally, however, from experiments the phenomenological spectral 
function Pph{s) is only known from threshold up to some energy Sq. Above this value, 
we shall use the theoretical expression pth{s) also for the right-hand side. In the case of 
the 5-mesons, we approximate the phenomenological spectral function by the pole of the 
lowest lying hadronic state plus the theoretical spectral function above the threshold sq, 

pph{s) = m%fl5{s - m]^) + e{s - SQ)pth{s) . (1.8) 

This is legitimate if sq is large enough so that perturbation theory is applicable. The 
central equation of our sum-rule analysis for Jb then takes the form: 

so 

mUl = J e(™^-^)/>,(s)ci.. (1.9) 



Besides the sum rule of eq. ( |1.9|) , in our numerical analysis we shall also utilise a second 
sum rule which arises from differentiating eq. ( [l.T] ) with respect to 1/u: 



d 



d{l/u) 



oo oo 

J se~''''p{s)ds = m|/|e-'"^/" + J s e-'''' pth{s) ds . (1.10) 



so 



Taking the ratio of the sum rules of eqs. ( |1.10| ) and (|1.9| ), the decay constant drops out. 



and, as far as the phenomenological side is concerned, we end up with a sum rule which 
only depends on the heavy meson mass uib- In our numerical analysis, this additional sum 
rule will be used to fix the continuum threshold Sq from the experimental value of m^. 
The resulting sq is then used in the fs sum rule of eq. (pTQ]). 

In section 2, we give the expressions for the perturbative pseudoscalar spectral function 
up to the next-next-to- leading order in the strong coupling, and in section 3, the non- 
perturbative condensate contributions are summarised. Section 4 contains our numerical 
analysis of the sum rules. Finally, in section 5, we compare our results to previous deter- 
minations of /b in the literature and we present an estimate of the hadronic 5-parameter 
in the 5-meson system Bb- 
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2 Perturbative spectral function 



In perturbation theory, the pseudoscalar spectral function has an expansion in powers of 
the strong couphng constant, 

p{s) = p(^\s)+p^'\s)a{fia)+P^'\s)a{f,af + ... , (2.1) 

with a = as/ 71- The leading order term p^^\s) results from a calculation of the bare 
quark- ant iquark loop and is given by 



pio) 



» = 0(M + m)^s(l-^) . (2.2) 



For the moment, we have only kept the small quark mass m in the global factor [M + m)^ 
and have set it to zero in the subleading contributions. Higher order corrections in m up 
to order will be discussed further below. 

Our expressions for the spectral function always implicitly contain a 6'-function which 
specifies the starting point of the cut in the correlator \l/(s). Although generally, we prefer 
to utilise the MS mass, in order to have a scale independent starting point of the cut, 
in this case we chose the pole mass Mpoie- Modulo higher order corrections, it is always 
possible to rewrite the mass in the logarithms which produce the cut in terms of the pole 
mass such that the 6'-function takes the form 9{s — M^^^J. 

The order correction for the two point function \E'(s) was for the first time correctly 
calculated in ref. [^], keeping complete analytical dependencies in both masses M and m. 
Further details on the calculation can also be found in ref. ||2^. From these results it is a 
simple matter to obtain the corresponding imaginary part: 

p(i)(s) = ^^CF(M + mfs{l-x)\(l-x)\4:L2ix) + 2lnx\n{l-x) (2.3) 
Idtt^ l_ L 

1 1 
- (5 - 2x) ln(l - x) J + (1 - 2a;)(3 -x)lnx + 3(1 - 3a;) In ^ + -(17 - 33a;) 

where x = M'^/s and L2{x) is the dilogarithmic function [^]. The explicit form of the 
first order correction is sensitive to the definition of the quark mass at the leading order. 
Eq. (|2.3|) corresponds to running quark masses in the MS scheme, M{prn) and m{prn), 
evaluated a the scale pm- 

The term proportional to In pf^/M^ cancels the scale dependence of the mass at the lead- 
ing order, reflecting the fact that p{s) is a physical quantity, i.e., independent of the renor- 
malisation scale and scheme. Transforming the quark mass into the pole mass scheme]^ the 
Explicit expressions for the relation between pole and MS mass are collected in Appendix B. 
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resulting expression becomes scale independent and of course agrees with eq. (4) of 0. As 
shall be discussed in more detail below, however, the perturbative corrections to Jb in the 
pole mass scheme turn out to be rather large and we refrain from performing a numerical 
analysis of the sum rule in this scheme. Therefore, our expressions for the spectral function 
will only be presented in the MS scheme. 

The three-loop, order correction p^'^\s) has only been calculated very recently by 



Chetyrkin and Steinhauser |Tq, 17] for the case of one heavy and one massless quark. 



A completely analytical computation of the second order two-point function is currently 
not feasible. However, one can construct a semi-numerical approximation for p^'^\s) by 
using Pade approximations together with conformal mappings into a suitable kinematical 



variable ||2^,^. The input used in this procedure is the knowledge of eight moments for 
the correlator for large momentum x ^ 0, seven moments for small momentum x —>■ oo, 
and partial information on the threshold behaviour x —>■ 1. In our analysis, we have made 
use of the program Rvs.m which, contains the required expressions for p^'^\s) and was kindly 
provided to the public by the authors of HlGilT^. 

In ref. the pseudoscalar spectral function p(s) has been calculated in the pole 

mass scheme. Thus we still have to add to p^'^\s) the contributions which result from 
rewriting the pole mass in terms of the MS mass. The two contributions Aip'-^-* and A2p^'^^ 
which arise from the leading and first order contributions, respectively, are given by 



(2)/ 



A2P 



(2)/ 



87r2 
87r2 



(M + my s (3 - 20x + 21x^) r^^ - 2(1 - x)(l - 3x) r^^ , (2.4) 



(2) 



Cf{M + mf sr^^U (1 -x)(l -3x) 4L2(x) + 21nxln(l 



x) 



- (1 - x){7 - 21a; + 8x^) ln(l - a;) + (3 - 22a; + 29x'^ - 8x^) Inx (2.5) 



i(l_a;)(15_31a;)| 



Explicit expressions for the coefficients rm and r™'* can be found in Appendix B. Further- 
more, in ref. |ll6|,|T^ the renormalisation scale of the coupling pa was set to Mpoie. Since in 
our numerical analysis we plan to vary the scale pa independently from pm, the contribu- 
tion which results from reexpressing a{M) in terms of a{pa) in the two-loop part needs to 
be included as well. 

Close to threshold, in the pole mass scheme, the pseudoscalar spectral function behaves 
as v'^{aslTiv)^ where v = {1 — x)/{l + x) at any order k in perturbation theory. This 
behaviour, however, does not persist in the MS scheme, where for each order, an additional 
factor of 1/v is obtained, such that the order correction goes like a constant for v — >■ 
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0. Nevertheless, as we will see in more detail below, numerically the corrections for the 
integrated spectral function show a much better convergence than in the pole mass scheme. 

Let us now come to a discussion of the corrections in the small mass m. At the leading 
order in the strong coupling and up to order m^, they can, for example, be found in ref. 



X Nr , f . X 9 (1 + x) Mm^ (1 — 2x — x'^) m , 

^\s) = -|(M+m)^ 2(l-x)Mm-2m^-2 ^ + \,_.,. (2.6) 



The somewhat bulky expressions for the first order as correction can be obtained by ex- 



panding the results of p2|,|23| in terms of m and have been relegated to Appendix C. 
Numerically, the size of the order «<, corrections increases with increasing order in the 
expansion in m. However, even for the case of Bg the mass corrections in become 
negligible before the perturbative expansion for these corrections breaks down. 

In the process of performing the expansion of the results of |2^,^ in terms of m, it is 



found that starting with order logarithmic terms of the form In m appear in the expan- 
sion. They are of infrared origin, and in the framework of the operator product expansion 
it should be possible to absorb them by a suitable definition into the higher dimensional 
operator corrections, the vacuum condensates. If the operator product expansion is per- 
formed in terms of non-normal ordered, minimally subtracted condensates rather then the 
more commonly used normal ordered ones, the mass logarithms indeed disappear 



3 Condensate contributions 

In the following, we summarise the contributions to the two-point function coming from 
higher dimensional operators which arise in the framework of the operator product expan- 
sion and parametrise the appearance of non-perturbative physics, if the energy approaches 
the confinement region. Here, we decided to present directly the integrated quantity u'^{u) 
because the spectral functions corresponding to the condensates contain (5-distribution con- 
tributions. 

The leading order expression for the dimension-three quark condensate is known since 
the first works on the pseudoscalar heavy-light system 



(M + mfM{qq) e"^'/" 



1 



(3.1) 



To estimate higher order mass corrections in our numerical analysis, we have included the 
corresponding expansion up to order m? . From the mass logarithms of the perturbative 



order as and correction, it is a straightforward matter to also deduce the first order 
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correction to the quark condensate since the mass logarithms must cancel once the quark 



condensate is expressed in terms of the non- normal ordered condensate |27-29]. We were 



not able to find the following result in the literature and assume that it is new: 

u^^{u) = ^C7.a(MWM(,-g){r(0,^)-[l+(l-^)(ln^ + l)]e--V"}, 

(3.2) 

where r{n,z) is the incomplete F-function. Again, the term In/x^/M^ cancels the scale 
dependence of the mass at the leading order. 

The next contribution in the operator product expansion is the dimension-four gluon 
condensate. Although its infiuence on the heavy-light sum rule turns out to be very small, 
we have nevertheless included it in the analysis. The corresponding expression for the 
Borel transformed correlator is given by 

„$W(„) = l(M + m)=^(aFF)e-^^'/". (3.3) 

In some earlier works on the pseudoscalar sum rule this contribution appears with a wrong 
sign although of course this has negligible infiuence on the numerical results. 

The last condensate contribution that we consider in this work is the dimension-five 
mixed quark gluon condensate which still has some infiuence on the sum rule since it is 
enhanced by the heavy quark mass. Again here the result is well known from the literature 
and we just cite it for the convenience of the reader: 

„$™ („) = - (M + raf M(i££l> (i _ Ml) ,-MV« . (3.4) 

We have checked explicitly that the contribution of the next-higher dimensional operator, 
the four-quark condensate, is extremely small, and thus have neglected all higher dimen- 
sional operators. The corresponding results for the condensate contributions to the sum 



rule of eq. (|1.10|) can be calculated straightforwardly by differentiating the above expres- 



sions with respect to 1/u. 



4 Numerical analysis 

In our numerical analysis of the pseudoscalar heavy-light sum rule, we shall mainly discuss 
the values of our input parameters, their errors, and the impact of those errors on the values 
of fs and fs^- To begin, however, let us investigate the behaviour of the perturbative 
expansion. 
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As was already mentioned above, in the pole mass scheme the first two order as and a1 
corrections to ^{u) are of similar size than the leading term, thus not showing any sign of 
convergence. For central values of our input parameters and a typical value m = 5 GeV^, 
the first order correction amounts to 78% and the second order to 85% of the leading term. 
To be consistent with the perturbative result for p(s), we have used m^"^*^ = 4.82 GeV, 
which results from relation (|B.4|) up to order a^. Because of the large corrections, we shall 
not pursue an analysis in the pole mass scheme any further. On the contrary, in the MS 
scheme for /i^ = /ia = f^^b and m = 5 GeV^, the first and second order corrections are 11% 
and 2% of the leading term respectively, while at /i^ ~ 4.5 GeV the second order term 
vanishes entirely. Hence, in the MS scheme the perturbative expansion converges rather 
well and is under good control. 

240 I — • — • — ■ — ■ — \ — ■ — ■ — ■ — ■ — \ — • — • — ■ — • — \ — • — • — • — • — ^ — • — • — • — ■ — I 



230 : 
220 - 




u [GeV'] 

Figure 1: fBd^sa function of the Borel parameter u for different sets of input parameters. 
Solid line: central values of table 1; long-dashed line: nih^mf,) = 4.16 GeV (upper line), 
mb^rrib) = 4.26 GeV (lower line); dashed line: fim = 3 GeV (lower line), = 6 GeV 
(upper line). 

In figs. 1 and 2, as the solid lines we display the leptonic decay constants fs and fs, 
for central values of all input parameters which have been collected in tables 1 and 2, 
as a function of the Borel variable u. For "U ^ 4 GeV^ the power corrections become 
comparable to the perturbative term, whereas for u > 6 GeV^ the continuum contribution 
gets as important as the phenomenological part. Thus a reliable sum rule analysis should 
be possible in the range roughly given by 4 GeV^ ^ ^ ^ 6 GeV^. In this region we extract 
our central results fs = 210 MeV and fs, = 244 MeV. 
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Figure 2: fs^ as a function of the Borel parameter u for different sets of input parameters. 
Solid line: central values of table 2; long-dashed line: mb(m{,) = 4.16 GeV (upper line), 
mbinih) = 4.26 GeV (lower line); dashed line: /i^ = 3 GeV (lower line), /im = 6 GeV 
(upper line). 



As an additional input parameter the continuum threshold Sq is required. This param- 
eter can be determined from the ratio of the sum rules of eqs. ( |1.10| ) and ( |1.9D , which only 
depends on the heavy meson mass. To this end, for a certain set of input parameters, Sq 
is tuned such as to reproduce the Particle Data Group values for and |RD[ in the 
stability region (a minimum in this case) of the ratio of sum rules. In tables 1 and 2, we 
also present the resulting values for sq and the corresponding location uq of the minimum 
of the iriB sum rule. For central values of all input parameters, we obtain sq = 33.6 GeV^ 
and Mo = 5.6 GeV^ for the _B-meson, as well as sq = 35.5 GeV^ and uq = 5.1 GeV^ for 
the Bg-meson. In fig. ^, we show the resulting tub and m^^ as a function of u for central 
input parameters. As can be seen from this figure, in the stability region, the sum rule 
reproduces the physical heavy meson masses which are indicated as horizontal lines. Our 
results for fs and fs^ are then extracted at Uq, around which also the sum rules for the 
decay constants are most stable and display an inflection point. 

The dominant source of uncertainty for the decay constants is the error on the bottom 



quark mass mi,. For this value we have taken an average over recent determinations |[3TH39 
which results in mb{mb) = 4.21 ±0.05 GeV. The error on rrib has been chosen such that all 
individual results are included within one standard deviation. The corresponding variations 
of and /^^ are displayed as the long-dashed lines in figs. 1 and 2, where the upper line 
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Parameter 


Value 


So [GeV^] 


uo [GeV^] 


A/b [MeV] 


mb{mb) 


4.21 ±0.05 GeV 


33.1 
34.2 


6.1 
5.2 


±15 




3.0 - 6.0 GeV 


33.5 
34.4 


6.8 
4.0 


±10 




3.0 - 6.0 GeV 


34.2 
33.1 


5.1 
6.2 


+2 
-1 


(gg)(2 GeV) 


- (267 ± 17 MeV)=^ 


33.9 
33.3 


5.7 
5.5 


±6 




2xO(a2) 
no 0(02) 






±2 




0.1185 ±0.020 






±1 


{aFF) 


0.024 ±0.012 GeV^ 






±1 


ml 


0.8 ± 0.2 GeV^ 






±1 



Table 1: Values for all input parameters, continuum thresholds sq, points of maximal 
stability Uq, and corresponding uncertainties for fs- 

corresponds to a lower value of m^ and the lower line to a larger m^. The impact of the 
variation of m^ on the error of fs and Jb^ has been quantified in tables 1 and 2. 

Another important source of uncertainty is the renormalisation scale /im- We have 
decided to vary fim in the range 3-6 GeV, with a central value fim = rrib. If /im is 
smaller than about 3 GeV, the perturbative corrections become too large and the expansion 
unreliable. As the dashed lines in figs. 1 and 2, we then show the corresponding results for 
fim = 3 GeV (lower line) and /i^ = 6 GeV (upper line). The uncertainties for fs and /^^ 
which result from /i^ are again listed in tables 1 and 2. To indicate the influence of even 
lower scales, let us briefly discuss the GeV. Here, we find sq = 38.4 GeV^ 

being rather large, as well as Uq = 2.9 GeV^ which is very small. At such a low uq, the 
perturbative and operator product expansions are not very reliable. Nevertheless, the value 
for fs extracted at Uq turns out surprisingly close to our central result, such that the error 
estimate of table 1 is more conservative. The variation of fia, on the other hand, only has 
a minor impact on the error of Jb and Jb^ and is also given in tables 1 and 2. 

The present uncertainties in the remaining QCD parameters a^, the strange quark 
mass mg and the condensate parameters have much less influence on the errors of Jb and 
/b^. Thus let us be more brief with the discussion of these quantities. The current value 



of as{Mz) by the Particle Data Group, as{Mz) = 0.1185 ± 0.020 has been used, 
whereas our choice for the strange mass ms(2 GeV) = 100 ±15 MeV is obtained from 
two very recent analyses of scalar and pseudoscalar QCD sum rules The resulting 



mg is compatible to the determination from hadronic r-decays, as well as lattice QCD 
results [^2|-0. Besides the variation of as{Mz), in order to estimate the influence of 
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Parameter 


Value 


So [GeV^l 


no [GeV^] 


AfB [MeVl 


mb{mb) 


4.21 ±0.05 GeV 


34.8 
36.4 


5.4 
4.8 


^16 




3.0 - 6.0 GeV 


35.2 
37.2 


6.2 
3.6 


+8 
_9 




3.0 - 6.0 GeV 


36.2 
34.9 


4.7 
5.5 


+1 


{ss)/{qq) 


0.8 ±0.3 


35.9 
35.2 


5.3 
4.7 


±8 


{qq){2 GeV) 


-(267± 17MeV)3 


35.7 
35.3 


5.2 
4.9 


+5 
-4 


m,(2 GeV) 


100 ± 15 MeV 






±2 




2xC'(a2) 
no ©(ai) 






±3 




0.1185 ±0.020 






±1 


(aFF) 


0.024 ±0.012 GeV^ 






±1 


ml 


0.8 ± 0.2 GeV^ 






±1 



Table 2: Values for all input parameters, continuum thresholds Sq, points of maximal 
stability uq, and corresponding uncertainties for /b^- 



higher order corrections, we have either removed or doubled the known 0{a1) correction. 
The resulting uncertainty for the decay constants, however, turns out to be small. 

Our value for the quark condensate has been extracted from the Gell-Mann-Oakes- 
Renner relation |55] with current values for the up- and down-quark masses |5l]. The 



ratio (ss) / (qq) has been chosen such as to include results from refs. [P,^6H49[|.n The mixed 
quark-gluon condensate is parametrised by {gsqcrFq) = ml{qq) with ml being determined 



in ref. [0], and finally, for the gluon condensate we take a generous range which includes 



previous values found in the literature. All uncertainties for and fs^ resulting from 
these parameters are also listed in tables 1 and 2. Where entries for sq and Uq are missing, 
we have used the values corresponding to central input parameters. 

Adding all errors for the various input parameters in quadrature, our final results for 
the B and Bg meson leptonic decay constants are: 

fs = 210 ±19 MeV and = 244 ±21 MeV. (4.1) 

In the next section, we shall compare these values with previous QCD sum rule and lattice 
QCD determinations. 

''We have not taken into account the very recent result (ss) / (qq) = 1.7 [Q, obtained in the framework 
of chiral perturbation theory (xPT), which would lead to fs^ — 270 MeV. In xPT, the value of the quark 
condensate depends on the subtraction procedure employed, and it is not clear how these results relate to 
(qq) in the MS scheme. The large value obtained in can almost be excluded on the basis of our Jb, 
together with independent lattice results for the ratio fsj fs (see below). 
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Figure 3: niB (solid line) and m^^ (dashed line) as a function of the Borel parameter u 
for central input parameters. The horizontal lines indicate the corresponding experimental 
values for these quantities. 

5 Conclusions 

The only truly non-perturbative method to compute hadronic matrix elements is QCD 
on a space-time lattice and thus it is very interesting to compare our findings to the 
corresponding results in lattice gauge theory. For the leptonic heavy meson decay constants, 
they have been compiled in a recent review article by Bernard [^.0 Taking into account 
dynamical sea quark effects and estimating the corresponding uncertainties, his world 
averages read: 



B 



200 ± 30 MeV and 



Ib 



1.16 ±0.04. 



(5.1) 



The lattice value for fs is in good agreement with our result of eq. ( [4 .11) , and also our 



ration /bs/ Ib = 1-16 turns out to be perfectly consistent with (|5.1| ). Nevertheless, due to 
sizable discretisation errors on the lattice, in our opinion, at present the QCD sum rule 
determination of the decay constants is more precise. 

We now come to a comparison with recent QCD sum rule results for fs and Jb,- The 
status of sum rule calculations of Jb in the pole mass scheme has been summarised in 
the review article |]14| with the result Jb = 180 ± 30 MeV. Although roughly 15% lower. 



within the errors this result is compatible with our value (|4.1| ). However, due to the large 



^See also ref. 
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perturbative corrections in the pole mass scheme, and the strong dependence on the bottom 
quark mass which in [14] was taken to be rri^°^^ = 4. 7 ±0.1 GeV, the theoretical error is not 
controlled reliably. Let us remark that the order correction in the relation between MS 



mass and pole mass alone gives a shift of m?°^^ by roughly 200 MeV [p^-|56[|. Typical results 



for fs^ turn out to be about 35 MeV higher than fs so that the difference between 
fs^ and fs is in agreement to our result. Our result for fs is also completely compatible 



with the very recent analysis of ref. ||20|, which was performed in the framework of HQET 
and resulted in Jb = 206 ± 20 MeV, suffering however from the problems of the pole mass 
discussed above. 

After submission of our work to the e-print archive, an independent analysis of the 



heavy-light meson sum rules by Narison [57| was published, which also employs the heavy 



quark mass in the MS scheme. For the convenience of the reader, even though ref. [5^ 
appeared later, we have been asked by our referee to nevertheless comment on this analysis. 
The main difference to our analysis lies in the fact that in ref. the bottom quark mass 
is extracted from the sum rule for rriB, with the result mbirnh) = 4.05 ± 0.06 GeV. We 
have checked that for this value of m?, one needs sq = 37.5 GeV^ to reproduce tub, and 
finds a stability region around uq = 4.3 GeV^. Inserting these parameters into the /b sum 
rule, we obtain Jb = 270 MeV, in conflict to our result ( |4.1| ). We are able to reproduce the 
value quoted by Narison, /b = 205 MeV, at m = 2.7 GeV^, which roughly corresponds to 
his preferred r = 1/m value. Around this u, however, the fs sum rule is unstable, casting 
doubts on the procedure of also demanding stability in the continuum threshold sq, besides 
the M-stability. Furthermore, the rather low value of rrib compared to our world average 
presented above, as well as the very high value of Sq, indicate that the pseudoscalar sum rule 
is not a good place to determine rrih. On the other hand, our investigation demonstrates 
that perfectly compatible results are obtained with a more standard value for nib. The ratio 



IbJIb has been calculated by the same author in ||5^ with the result /bs/ fs = 1.16±0.05, 

in agreement to our findings. 

The heavy-meson decay constant Jb plays an important role in the mixing of neutral 
and mesons. The relevant hadronic matrix element can be expressed as |59 



{B'\Qab=2\B') = ^BBflml, (5.2) 

where Qab=2 is the scale invariant four-quark operator which mediates B^-B^ mixing and 
Bb is the corresponding scale invariant 5-parameter which parametrises the deviation of 
the matrix element from the factorisation approximation. In the factorisation approxima- 
tion, by definition we would have Bb = 1- The combination ^/BbJb can be extracted from 
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an analysis of experimental data on B^-B^ mixing together with additional inputs which 
determine the matrix elements of the quark mixing or Cabibbo-Kobayashi-Maskawa ma- 
trix. A very recent analysis then yields ^/BbIb = 236 ± 35 MeV |^|^. Taking together 



our result for fs and the quoted value for ^/Bb/b, we are in a position to give an estimate 
of the scale invariant 5-parameter Bb, which reads 

Bb = 1.26 ±0.45. (5.3) 

For simplicity we have assumed Gaussian errors in both input quantities. The result again 
is in very good agreement to corresponding determinations of Bb on the lattice which gave 
Bb = 1-30 ± 0.12 ± 0.13 although our error in this case is bigger. 

To conclude, in this work we have presented a QCD sum rule determination of the 
leptonic heavy- meson decay constants /b and J'b^- Due to large perturbative higher order 
corrections, an analysis in terms of the bottom quark pole mass appeared unreliable. On 
the contrary, employing the heavy quark mass in the MS scheme, up to order the per- 
turbative expansion displays good convergence and a reliable determination of and /b, 
turned out possible. Our central results have been presented in eq. (|4.1| ), where the dom- 
inant uncertainty arose from the present error in the bottom quark mass mf,(mb). Taking 
into account independent information on y/Bsfs from B^-B^ mixing, we were also in a 



position to give an estimate on the 5- meson S-parameter Bb in eq. ( |5.3|) . All our results 
are in very good agreement to lattice QCD determinations of the same quantities. Further 
improvements of our results will only be possible if the dominant theoretical uncertainties 
could be reduced. This would require a more precise value of the bottom mass, and a 
reduction of the renormalisation scale dependence, requiring the next perturbative order 
correction, which at present seems to be out of reach. 
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Appendices 



A The Borel transform 

The Borel operator Bu is defined by (s = — 



Bu = lim 

s,n— >oo 

s/n=u 



{n - 1)! {dsy 



(A.l) 



The Borel transformation is an inverse Laplace transform [^]. If we set 



f{u) = Bu f{s) , then f{s) 



u 



f{u) e-'/'^du , 



(A.2) 



In this work we just need the following Borel transform: 



a, 



1 



[X + s) 



1 



-x/u 



(A.3) 



Cases in which logarithms appear can be treated by first evaluating the spectral function 
and then calculating the dispersion integral of eq. (pTTp. 



B Renormalisation group functions 

For the definition of the renormalisation group functions we follow the notation of Pascual 
and Tarrach |6^, except that we define the /3-function such that f3i is positive. The 
expansions of /3(a) and 7(a) take the form: 



/3{a) 



and 7(a) = 7ia + 72^^ + 73^^ + • • • , (B.l) 



with 



Pi 



and 



UCa - 4:Tn 



71 



P2 



72 



1 

12 

Cf 
48 L 



17Cl - IOCaTuj - GCpTrij 



97Ca + 9Cf - 20Tn 



f 



The relation between pole and running MS mass is given by 



M(/i, 



pole 



1 + a(/i„) r(^) (flrn) + a{Haf (/^a, flm) + ■■■ 



(B.2) 
(B.3) 

(B.4) 
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where 



,(1) 



,(2) 



(1) , /^m 

^m,o - 71 In 



.(2) 
m,0 



M(/im) 
72 + (71 - Pi) r% 



(B.5) 



, A'-m I 7l / a \ ^ 2 f^iri 



7i + /^i In — 



,(1) 



(B.6) 



The coefficients of the logarithms can be calculated from the renormalisation group |^| 
and the constant coefficients r|^ o ^-nd r^^Q are found to be [|6^ ,|65 



.(1) 

m.O 



.(2) 
m.O 



(B.7) 



+ ^^^^(-^ + k(2) + |C(3) - ^C(2) ln2) + CpT(^ - ^C(2) 



384 2 



V4 2 



(B.8) 



C Mass corrections at order a. 



Below, we present the order mass corrections to the pseudoscalar spectral function which 
arise from expanding the results by [^,^ up to order m^, after the higher dimensional 
operators have been expressed in terms of non-normal ordered condensates: 



Cf{M + mf Mm {{l-x) AL2{x) + 21nxln(l - a; 



(C.l) 



2(4 - x) ln(l - x) + 2(3 - 5x + x^) Inx + 3(2 - 3x) In + 2(7 - 9a;) 

M2 



pS(-) 



Cf {M + mYm^ { - x) 4L2(x) +21na;ln(l - x 



(2 + x)(4-x)ln(l -x) + (6 + 2x-x^)lnx + 61n-^ + (8 -3x) }, (C.2) 



p2(^) 



^ .oMm^ \ , . ^, , ,^ , (9 + 8x-9x2) 
-^Cir(M + m)2 <^ 4^2(0;) +21na;ln(l -x) + ^ - 



-2 



87r 

(7 + 7x - 2x2 
(1 -a;) 



(1-x)^ 

, ^ ,(6 + 7x-2x2), ,(2-x2), 

In 1 - x) + 2^ ^ In X + 6)^ In ^ 

^ ' {^-x) (l-x)2 M2 



(C.3) 
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N 

^ CiT (M + m)' — <j 2L2(a;)+lna;ln(l 



X) 



(13-24x-27x2 + 2x=^), , (12-22x-27x2 + 2a;3) , 

m(l — xj + — m X 



+ 3 



2(1 

(4 - 12,r + ,r2 + 3x^) 
2(1- x)3 



In 



M2 



2(1 

(6 - (31,r + 15,r- + 11 x''^) 
4(1 



(C.4) 
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